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We study the low-energy bandstructure of armchair and small-bandgap semiconducting carbon 
nanotubes with proximity-induced superconducting pairing when a spiral magnetic field creates 
strong effective spin-orbit interactions from the Zeeman term and a periodic potential from the 
orbital part. We find that gapless Dirac fermions can be generated by variation of a single parameter. 
For a small-bandgap semiconducting tube with the field in the same plane, a non-degenerate zero 
mode at momentum fe = can be induced, allowing for the generation of topologically protected 
Majorana fermion end states. 

PACS numbers; 73.63.Fg, 74.45.+C, 74.70.Wz 



I. INTRODUCTION 



The electronic properties of single-wall carbon nan- 
otubes (CNTs) have been studied for almost two decades 
by now and are generally thought to be well understood.^ 
The electronic structure of a CNT on energy scales be- 
low hvp/R, with radius R and Fermi velocity vp — 
8 X 10''^ m/s, is captured by the low-energy approach, 
where one retains only the lowest transverse momen- 
tum bands and disregards trigonal warping. We fo- 
cus on small-bandgap semiconducting CNTs without 
primary gap, where the curvature-induced bandgap is 
Eg ex cos(36')/i?^ and typically of order meV in ex- 
periments. Here, 9 is the chiral angle, ^ and only 
armchair tubes, cos{39) = 0, are metallic. Two in- 
teresting developments concern spin-orbit interactions 
(SOI) and proximity-induced pairing correlations. First, 
Coulomb blockade spectroscopy experiments for ultra- 
clean CNTg2EI have confirmed the existence of the theo- 
retically expected but rather small curvature-induced in- 
trinsic SOI.BiEi Although a proposal exists to design tun- 
able SOI in graphene by deposition of suitable adatoms,^! 
this idea does not readily apply to CNTs. Second, 
proximity-induced superconductivity was experimentally 
demonstrated and has been usefuUj^xploited.l^tEl How- 
ever, only few theoretical studies^^^ have addressed the 
corresponding pairing terms in the CNT low-energy the- 
ory. 

In this work, we demonstrate that the combined ef- 
fects of strong SOI and proximity-induced superconduc- 
tivity in CNTs are responsible for emergent gapless Dirac 
fermions and Majorana bound states. The schematic set- 
up considered here is shown in Fig. [l] Apart from the 
spiral magnetic field used to generate strong spin-orbit 
couplings, the set-up is s imilar to recent proposals for 
semiconducting nanowires-'^^'^ In both cases, despite of 
the presence of the superconducting film generating the 
proximity-induced pairing in the wire (or CNT), one can 




Figure 1: (Color online) Schematic set-up of a carbon 
nanotube proximity-coupled to a superconductor, in a spa- 
tially dependent magnetic field with angles ip{y) and i?(j/), 
see Eq. ([2|. The chemical potential on the tube is assumed 
to be tunable via a nearby gate electrode. 



change the chemical potential /x in the wire via a gate 
voltage. Our first step below is to show that experimen- 
tally available^^ spiral magnetic fields offer strong effec- 
tive SOI in CNTs, see also Ref. 19. As has been discussed 
in Ref. HU] it is possible to apply a spiral magnetic field to 
the CNT despite of the presence of the superconductor. 
We then proceed with a symmetry analysis of all possible 
proximity-induced pairing terms in the low-energy theory 
for a CNT in contact to an s-wave BCS superconductor. 
Employing experimentally realistic parameters, we find 
generic gap closings upon variation of a single control 
parameter, usually corresponding to zero-energy states 
("zero modes") with finite quasi-momentum k. The gap 
closing can be probed experimentally by tunneling spec- 
troscopy. For small-bandgap semiconducting CNTs with 
magnetic field in the same plane, however, the zero mode 
is non-degenerate and at k = 0. Th is implies that Ma- 
jorana bound states (MBSsJAfiE^MHH form at the tube 
ends. Majorana fermions may be useful for topological 
quantum computation,'^ and their realization is now ac- 



tively pursued in many different material systems. Us- 
ing InSb nanowires, a set-up similar to Fig.^ but with in- 
trinsically strong SOI instead of the spiral magnetic field, 
has been stu died experimentally, and clear evidence for 
the predictecP^'i^ MBSs at the ends of the nanowire was 
reported from tunnel spectroscopyJ^Sl For quantum com- 
putation applications of MBS states,'2S! CNTs could yield 
an attractive alternative to InSb nanowires. Apart from 
the unique electronic and mechanical properties of CNTs 
and their wide availability, MBS networks required for 
braiding operations^ZISll could be implemented by cross- 
ing CNTs. Such CNT crossin gs hav e been experimentally 
realized already a decade ago.f^SlMl 

MBSs in CNTs have also been proposed in two other 
papers, relying on either the intrinsic curvature-induced 
SOl^"* or^^ on an electric-field induced SOI.^^ Our pa- 
per is more closely related to the Sau-Tewari proposaljl^ 
but differs in two important regards. First, the spi- 
ral magnetic field induces a much stronger SOI than 
the curvature-induced SOI in CNTs. Second, this field 
also generates an effective periodic potential along the 
CNT, which automatically breaks valley {K, K') degen- 
eracy and implies a greatly enhanced robustness of the 
zero mode sector from which we explicitly construct the 
MBS. As a consequence, while the valley mixing param- 
eter IS.KK' is essential in the Sau-Tewari proposal,"^ it 
plays no significant role in the parameter space relevant 
for MBS generation in our proposal. For completeness, 
we nonetheless keep I^kk' in our model Hamiltonian. 

The structure of this article is as follows. In Sec. |TT] 
we introduce the low-energy model for a CNT in a spiral 
magnetic field. We then continue in Sec. |III| with a gen- 
eral analysis of the proximity-induced pairing terms ap- 
pearing in the Hamiltonian. Results for the band struc- 
ture are presented in Sec. |IV| both for armchair tubes 
(Sec. IV A I and for small-bandgap semiconducting tubes 



IV B). In the latter case, in Sec. IV C we explicitly 



(Sec. 

construct the MBS wavefunction when the spiral mag- 
netic field is in the same plane as the CNT. The role 
of electron-electron interactions is briefiy addressed in 
Sec. |IVD[ and we conclude in Sec. |V] 



II. CNT IN A SPIRAL MAGNETIC FIELD 

In the absence of a superconducting substrate, the 
single-particle Hamiltonian for a straight CNT along the 
j/-axis readgi' 

Hq = -ihvpdyay + Egax + {evFR/2)ByrizCrx (1) 

9e^J■B^ 



+ ^KK'Vx 



■B(y)-s 



with Pauli matrices cr^-yz (fjx.y.z) in sublattice (valley) 
space, where the two sublattices correspond to the two 
carbon atoms forming the basis of the honeycomb lattice 
and the two valleys denote the K, K' points in the first 
Brillouin zone; Pauli matrices s 



functions near the K,K' points, i.e., a state with quasi- 
momentum fc = (with —ihdy — > k) sits right at the K 
(or K') point. In Eq. (fll we omit the intrinsic SOP^^ 
since the spatially dependent magnetic field B(2/) will 
generate much larger couplings. In the Zeeman term, 
the Bohr magneton is /is and we use g^ ~ 2 for the 
Lande factor. Note that the orbital field along the CNT 
(j/-)axis favors valley polarization in z-direction. We also 
added the standard valley mixing term ,121221 /S.kk' , which 
arises due to residual elastic disorder and favors valley 
polarization in x-direction. 

Writing the magnetic field in polar coordinates, 
cf. Fig.g 



B = B{y) 



cos[Lp{y)] sin[i9(y)] \ 
sin[(p(j/)] sin[i9(2/)] 
cos[d{y)] / 



we next perform a unitary transformation ,12E01 



u{v) 



62 



'd[s\n{ip)sx---cos{ip)sy] 



(2) 



(3) 



aligning the local spin quantization axis with the mag- 
netic field direction. The unitarily transformed Hamil- 
tonian, Hq = W HqU with Eq. (flj, then contains the 
effective SOpl 



hvf 



(4) 




act m spin space. 
This Hamiltonian acts on slowly varying Bloch envelope 



x,y,z 



H' = —ihvpayll'dyU = ——<jySL[y) ■ s 



cos Lp sin "d 
sin ip sin d 
1 — cos d 



Here we consider spiral magnetic field configurations, 
where ^{y) = y/\ in Eq. (pi and both the field strength 
B and the angle Lp are constant. The case </? = (field 
in a plane perpendicular to CNT) could, for instance, 
be realized using the hyperfine nuclear fields discussed 
by Braunecker et al.^^ and the resulting SOI formally 
coincides with the electric-field induced SOI studied in 
Ref. 131] where we also recover their helical state solu- 
tions. The case <p = 7r/2 (field in the same plane as 
the CNT) has been realized experimentally using mag- 
netic superlattices.13 Note that for arbitrary (yS, the en- 
ergy scale 6 — hvp/{2X) associated to the "pitch" length 
A sets the effective SOI strength. Using a typical value 
A « 250 nm, 6 is several orders of magnitude larger than 
the previously discussed spin-orbit couplings in CNTs. 

When including proximity-induced pairing, it is 
convenient to work with Nambu spinors, ^'{y) ~ 
(V'|,'0|,V';j-V't); where H = /dy*^'H*/2. In a spi- 
ral magnetic field, using Eq. ^ and Pauli matrices r^^y^z 
in particle-hole space, the CNT Hamiltonian reads 

y. = [-ihvpcrydy + Egax + Asm{y/X)rizax - fj] r^ 
- fisBsz + 5ay[sm{ip)s,r: - cos{ip)sy]Tz 

+ ^KK'VxTz +'Ha, (5) 

where we added the chemical potential /i and a 
proximity-induced pairing term "Ha- Note that an or- 
bital field in tube direction causes a periodic potential 



with amplitude 



A 



evp 



RB sin (p. 



(6) 



Experimentally, it turns out that this estimate for A in 
some devices is enhanced by a factor 2 to 3 of unknown 
origin.HS Without SOI and superconducting pairing, pe- 
riodic p otenti als have recently been studied theoretically 
in CNT^24|35| ^^^ graphene.lSfil 

At this point, we briefly comment on the time-reversal 
symmetry (TRS) properties of "H. TRS requires that H 
commutes with the anti-unitary time-reversal operator 
T = iSyTjxC, where C denotes complex conjugation and 
T^ = —1. Requiring TRS for Ha, the only terms in 
Eq. (Is I violating TRS are the Zeeman term oc B and the 
orbital flux ex A. 



III. PROXIMITY-INDUCED 
SUPERCONDUCTIVITY 

Next we discuss the proximity effect due to an s-wave 
singlet superconducting substrate. To that end, we first 
write down all s-wave singlet pairing terms in the CNT 
which are consistent with TRS. With unity operator tjq 
(tJo) in valley (sublattice) space, we obtain the "zoology" 
of allowed pairing terms. 



Ha 



i—0,x 



Ai{cosxi Vx +sinxj r]o)a^Tx 



(7) 



-f- A2(cosx2 Vx + sinx2 r]o)<^yTy 

+ A3(cosx3 cTx + sinx3 cro)VyTx + /^i-qyCJyTy. 

Terms ex a^ (oc rj^) correspond to different substrate cou- 
plings for the two sublattices (valleys) and have been 
omitted in Eq. (l7|. In a generic situation, such asym- 
metries are extremely smalP^ except for very thin CNTs, 
where the low-energy approach does not apply in any 
case. We also did not include TRS-invariant terms 
oc Sx,y,z ill Eq. (Q, e.g., Sx<JyTx or SzTy. Such terms 
describe triplet pairing in the CNT, which cannot be gen- 
erated by coupling to an s-wave superconductor. 

The large number of parameters in Eq. (Q can now 
be greatly reduced by resorting to physical arguments. 
The simplest and most likely dominant contribution to 
the proximity effect comes from intra-sublattice (same 
orbital) pairing: In a microscopic lattice model. Ha — 
^j Ci^-|-c^,iAi, where i runs over all atoms in the hon- 
eycomb lattice. This pairing mechanism corresponds to 
the three terms oc a^ in Eq. ([7]). Another possibility 
comes from pairing between nearest-neighbor atoms (dif- 
ferent sublattices). This mechanism explains all remain- 
ing terms [oc <Jx,y\ in Eq- y\, but it has a much smaller 
amplitude because the two sublattices are not commen- 
surate. This implies that the 2fci?-oscillatory anomalous 
Green's function in the superconductor basically aver- 
ages out. Neglecting these subleading contributions, we 
are left with the terms oc Tx,rixTx and 7]yTx in Eq. (Q. 



Next, note that the r]xTx and rjyTx terms, which connect 
different valleys, are unitarily equivalent in the absence 
of KK' mixing. Since the gap closings reported below 
are also found for Akk' = 0, we omit, say, the r]yTx term 
in Ha- We then arrive at 



Ha == A(cosx nx + sinx ?7o)t:e, 



(8) 



with the proximity-induced gap A and the "pairing an- 
gle" X- This angle interpolates between pure inter-valley 
(x — 0) and intra-valley (x — t^ /2) pairing. The actual 
value for x depends on how strongly rotational symme- 
try around the CNT axis is broken by the presence of 
the substrate. If rotational symmetry stays intact, dif- 
ferent valley s form time-reversed partner states and have 
to be paired, I^EZI resulting in x = 0. On the other hand, 
TRS-invariant intra-valley pairing (x = ^e/^) dominates 
for strongly broken rotational symmetry.131 

Let us now briefly consider the case without the spiral 
magnetic fleld, A = B = 5 = G in. Eq. (l5]), where rjx — 
?7 = ± is conserved. Introducing the two proximity gap 
scales 



A, 



|smx + ?7Cosx| 



(9) 



the dispersion relation follows from 



El±{k) 



7]Akk' - M ± \l{hvFky + El 



At 



(10) 

When only one pairing term is present (x = or x = 
7r/2), the two gap scales coincide, A± = A, and the two 
pairing mechanisms cannot be distinguished unless the 
magnetic fleld is also present. Otherwise, however, the 
pairing angle is detectable since then A+ ^ A_. In fact, 
the dispersion relation may become gapless for rj — — 
and X = 7r/4. 



IV. RESULTS 

Analytical diagonalization of the Hamiltonian (pi with 
Ha in Eq. (|8| is not possible except for special cases, and 
in general we have to resort to numerics. Bloch's theorem 
impHes that eigenstates for the nth energy band, En{k), 
are of the form 



*fc,«(2;) 



— gifey 






mSZ f 



'i'V/>^^i 



) 



(11) 



with the multi-index v = (a, rj., s, r), where (Tz|(7 = ±) = 
cr|cr) (and so on). The quasi-momentum k is taken in the 

(n) 

flrst Brillouin zone, —1 < 2k X < 1, and ^rri'iy determines 
the normalized eigenstate. In this basis, all matrix ele- 
ments of H in Eq. (jsl except for the periodic potential 
oc A are diagonal in m. The A term couples m and to ± 1 
states, and diagonalization of the resulting Hamiltonian 
matrix yields En{k) and the eigenstates ( 11 ). We always 
flnd En{~k) — En{k) and thus show only half of the 
Brillouin zone below. 



The model parameters in Eqs. (|5| and (|8| are chosen as 
follows. We take an effective SOI scale 5 = 2 meV corre- 
sponding to magnetic pitch length A « 250 nm. This is a 
typical value for magnetic super lattices, see Ref. [TH) The 
Zeeman scale is ^bB = 0.5 meV (for _B « 5 T), and tak- 
ing into account Ref. |32j the orbital field implies the peri- 
odic potential amplitude A = Ai^sunp, see Eq. ([6|, with 
Aq = 2 meV (for i? sa 1 nm). For the proximity-induced 
gap, A = 0.3 meV is appropriate for Nb substrates^ 
The KK' mixing scale is taken as /S.kk' = 0.2 nieV,'2'22l 
which is a typical order-of-magnitude value for this phe- 
nomenological parameter. We note again that in con- 
trast to the proposal of Ref. O /S.kk' is not necessary 
for the MBS generation reported in Sec. |IVD[ Next, for 
the bandgap we choose either Eg = (armchair case) or 
Eg = 3 meV (small-bandgap semiconducting CNT). For 
the field angle (p, we take mostly ip = or (p = tt/2, but 
we have checked that small deviations from these values 
do not cause qualitative changes (see also below for a 
more detailed discussion). This leaves us essentially with 
two free parameters, namely the chemical potential fj, and 
the pairing angle x- (The pairing angle, however, is hard 
to change in an actual experiment.) 



A. Armchair CNT^ 

We begin with the armchair case. Eg = 0, where 
Vx = V — ^ '^^ conserved. For field angle ip = 0, the 
spectrum for Eq. (Is]) with (Is]) can be obtained analyti- 
cally by squaring H twice. Noting that also ay = a ~ ± 
is conserved. 



E: 



ri,a,± 



(k) 



k7](T 



U^bB) 



At 



(12) 



± 2 



k'qcr 



[<52 



{^iBBY] + {^1BB\Y 



with Ckrja = CTVph — ^ + riAKK' and A^ in Eq. ([9]). Equa- 
tion ( 12 1 does not permit zero-energy solutions for finite 



S and A^, and there is no gap closing for this configura- 
tion. 

Turning to p = 7r/2, numerical solution at chemical 
potential /i = yields the energy bands in Fig. [2j No- 
tably, the gap closes for x — 1-01, see inset of Fig. |2] 
with two pairs of zero modes at finite momenta k. Note 
that for both values of k, each zero mode is still twofold 
degenerate because of its partner state at — /c. The linear 
(massless Dirac fermion) dispersion is clearly visible in 
Fig. [2] Allowing also for a finite but small value of /i, 
we find that only one of the gap closings seen in Fig. [2] 
persists. 



B. Small-bandgap semiconducting CNTs 



Next we study the case of a small bandgap; for con- 
creteness. Eg = 3 meV. For field angle (^ = 0, we find just 
a single zero mode at finite k, along with its partner state 
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Figure 2: (Color online) 
case with field angle ip = 
pairing angle x ~ 0.327r 
positive k where the gap 
states with E„{—k) = E„ 
first Brillouin zone. Inset: 
diamonds give numerical 
the eye only. 
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Dispersion relation for the armchair 
it/2, chemical potential /i = 0, and 

= 1.01. Note that there are two 
vanishes, plus the respective A; < 
(fc) , where we show only half of the 

Gap closing as a function of x- Red 
results, the dotted line is a guide to 




Figure 3: (Color online) Low-energy dispersion relation for 
a CNT with bandgap Eg = 3 meV, field angle p = n/2, 
chemical potential fj, = 2.86 meV, and pairing angle x = 0.14. 
Note the linear dispersion near fc = 0, where the velocity is 
V ~ 0.15uf. The inset shows the dispersion also for higher 
energies. 



at —k. To give concrete numbers, a finite-fc zero mode 
was found for /z = 2.96 meV and x — 0.76. This behavior 
is qualitatively similar to the armchair case with p = 7r/2 
and finite /i. 

For a small-bandgap CNT and field angle p = 7r/2, 
we encounter a remarkably different situation with only 
a single non-degenerate zero mode at /c = 0, see Fig. [3] 
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Figure 4: (Color online) Band gap AEq (at fc = 0) for the 
parameters in Fig. [S] but with varying field angle (p. Red 
diamonds indicate numerical results, the dotted curve is a 
guide to the eyes only. 



Hr, 



-ihvTzdy + ayfx 



(13) 



The second term is due to the spatial variation of the 
chemical potential, where we find r^ ~ fx in the zero- 
mode basis by using the numerically obtained eigenvec- 
tors <&(='=). Writing $ = c+$^+) -|- c_$(~) (where c± are 
complex numbers), the state 



*(?/) ex e-"a'/2'^"$, fy<^ = -$, 



(14) 



then yields a zero-energy Majorana fermion solution, 
"Hp^ = 0, localized near y = 0. This state is topolog- 
ically protected by the gap to the next excited state. For 
the parameters of Fig. |3] this ensures MBS robustness 
for temperatures T < 1 K. 



D. Electron-electron interactions 



Only in such a non-degenerate case, single MBS forma- 
tion is possible. ^'^ Note that TRS has been broken by 
the applied magnetic field here; otherwise the M BS must 
have an overlapping time-reversed partner .ESBHl This gap- 
less state can be reached upon variation of x through a 
"magic angle," cf. inset of Fig.|2] An experimentally eas- 
ier route is to change, for fixed x- the chemical potential 
to the critical value through variation of a gate voltage. 
We stress that this zero-energy state can be reached by 
tuning a single parameter. 

We have also checked that the fc = gap closing ob- 
served in Fig. [3] stays robust against small changes of (p 
(or other model parameters). To illustrate this point, we 
show the fc = gap AEq as a function of (p in Fig. |4] 
Clearly, for ip near tt/2, we find AEq = 0, but sufficiently 
large deviations will destroy the zero mode. 



C. MBS construction 

For small-bandgap semiconducting CNTs with (p = 
it/2, the fc = zero mode in Fig. Is] and the gap closing 
and reopening as /i is varied through its critical value fic 
suggest that a MBy^lHSS] exists at the interface of regions 
with ^ < fj,c and ^ > fj,c. In practice, for fj, < iic, a MBS 
should then form at the CNT ends. To model this situ- 
ation, let us consider /i(j/) = fic ~ cty with a > 0, where 
a MBS is expected near y = 0. We explicitly construct 
the MBS wavefunction by first projecting the full Hamil- 
tonian (Isl to the Hilbert space spanned by the massless 



■H* 



(±) _ 



zblifc^ 



(±) 



ith 



O.Wvf 



,i±)i 



and ^f. (y) known numerically; the eigen- 



Dirac fermions 
from Fig. 3 

vectors ^m.l in Eq. (Ill are evaluated at fc = 0. With 
Pauli matrices f. 
by <&(+) and $' 



x,y.z acting in the Hilbert space spanned 
'', the projected low-energy Hamiltonian 



So far we have ignored electron-electron interaction 
effects beyond mean-field theory. Following the well- 
known fact that interactions destabilize the Fermi liq- 
uid phase in one dimension, it has been suggested ,33121 
and subsequently observed,B2Hl5l that CNTs display Lut- 
tinger liquid behavior. The Luttinger liquid phase is a 
strongly correlated phase, and therefore one should be 
careful in applying noninteracting theories to CNTs as 
done here. Importantly, the experiments in Refs. I43fl45l 
were performed without close-by metallic gates such that 
the long-range character of the Coulomb interaction was 
important. However, scanning tunneling spectroscopy for 
CNTs deposited directly on a metallic substrate did not 
show pronounced interaction phenomena (for a review, 
see Ref. 46), presumably due to the strong screening of 
the Coulomb potential by the substrate. Similarly, we ex- 
pect that the presence of the superconducting substrate, 
see Fig. [l] drastically reduces the effective interaction 
strength. The remaining weak interactions, on the other 
hand, are then not expected to destroy the MBS state 
or the emerging Dirac fermions discussed above. This 
topological stability against weak interactions has been 
discussed in detail in several recent works, "^^^-^ and we 
refrain from repeating their analysis here. 



V. CONCLUDING REMARKS 

In this work, we have studied the low-energy band- 
structure of CNTs with effective SOI (due to a spiral 
magnetic field) and proximity-induced pairing. As dis- 
cussed in Sec. |IVD| in the presence of the superconduct- 
ing substrate, interactions axe screened and our single- 
particle approach should be useful. Despite of the com- 
bined presence of a bandgap, the strong SOI, the orbital 
periodic potential, and the proximity gap, we find zero 
modes with associated massless Dirac fermions. The gap 



closings should guide future tunneling spectroscopy ex- 
periments: For parameters near a gap closing condition, 
topological end states will appear and can be observed 
as peaks in the dl/dV curve!^ For finite-fc zero modes, 
this corresponds to degenerate Majorana states at each 
end, while for the fc = gap closing, a single localized 
Majorana mode results. 
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